General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



Senl-Annual Progress Report 


on 

SUPERSONIC CO ELECTRIC-DISCHARGE LASERS 


NASA Grant NSG-2050 


fN 

(N 

9 


I 

vO 

P* 

SE 


'J1 

•9 r- 

o *- 

a <N 
s o 


for the period 

May 1, 1975 - October 31. 1975 


sO 

ry 

\ 

m 


Submitted to the 


National Aeronautics jnd Space Administration 

M 

Washington, D. C. 


u 

r-l Ift U 


a O' 
c «- 
s 

<tJ • 

O 

O ■ O 

a> o • 
u ^ 

M r- ^ 

as PI 

0(0 U 
r /1 ca I =2 
(B M 

gu Uj ,»«i 

(h •* >9 
^ aa as 

in 

w 

O 

^ os k ^ 
o> •< v • 

00 * u > 

vf) O O *H 

tn cu c 

9 M V s 
<- O OS 

1 t 

X U (A U 

0 *-( w o 

1 CB 0) tw 

< H w £S 
to U O' «J 
|i4 O 4J 
as •-) M CO 
W 04 — ' 


by 

Professor R> K* Hanson 
Professor M. Mitchner 
Mr. Alan Stanton 


High Temperature Gasdynamics Laboratory 
Department of Mechanical Engineering 
Stanford University 
Stanford. California 94305 



TABLE OF CONTENTS 

Page 

1.0 Sumnary 1 

2.0 Laser Modeling 2 

2.1 Probabilities for V-V Exchange in CO-N 2 and 

^ 2~^2 * ^ 

2.1.1 CO-N^ V-V Probabilities 2 

2.1.2 N 2 -N 2 Probabilities 6 

2.2 Calculations of CO/N 2 leaser Kinetics using C0ED2 ... 9 

2.2.1 Case 1; K 2 iCO = 90:10 12 

2.2.2 Case 2, N2.-C0 * 98:2 17 

3.0 Laborat '^ry Studies of Double Discharge 21 

4.0 Research Plan for November 1, 1975, to April 30, 1976 .... 23 

5.0 References 25 



I 


LIST OF FIGURES 


Figure 

1 RooT«-temperature probabilities for the CO-M 2 
V-V process CO(v) + NjCv-O) CO(v-l) + N 2 (v*l) . 

The calculated probabilities are fit to a value of 
1.9 X 10"^ for . 

oil 

2 Comparison of the calculated probabilities with the 

low-temperature experimental results of Stephenson 
and Mosburg for the V-V process C0(v=l) + N2(v«0) -► 
C0(v«0) + N2 (v*1) 

3 Calculated probabilities for the N 2 ~N 2 V-V exchange 

N2(v) + N2(v=0) - N2(v- 1) + N2(v*l) 

4 Normalized CO vibrational populations for Case 1, at 
four locations from the beginning of the discharge . . 

5 Maximum line-center gain coefficients for the CO laser 

transitions V_ + 1 V , Case 1 

L L 

6 P-branch and R-branch gain coefficients at line 
center for the v = 4 -»• v = 3 vibrational band. The 
distance from the beginning of the discharge is . 

5.5 cm 

7 Normalized CO vibrational populations for Case 2 . . . 

8 Maximum line-center gain coefficients for the CO 

laser transitions V + 1 -*■ V- , Case 2 . 

L L 


Page 

7 

8 

10 

13 

14 

16 

18 


19 


1.0 SUMMARY 


This report describes work carried out on NASA Grant NSG-20S0 during 
the period May 1975, to October 31, 1975. Our laser modeling activity 
has Involved addition of an option to the current program, C0ED 2, allow- 
ing N 2 as a second diatomic gas. This option is now operational and a 
few test cases involving N 2 /CO mixtures have been run. Results from these 
initial test cases are summarized in this reports 

In the laboratory, a new CW double-discharge test facility has been 
constructed and tested. New features Include: water-cooled removable 

electrodes, 0-ring construction to facilitate cleaning and design modifi- 
cations, increased discharge length, and addition of a post-discharge 
observation section. Preliminary tests with this new facility using N 2 
have already yielded higher power loadings than obtained in our first- 
generation facility. At present, power loading is limited by the power 

supply available, but a larger supply is on order and will be installed 
% 

shortly. Another test- sect ion modification, recently made and as yet un- 
tested, will permit injection of secondary gases into the cathode boundary 
layer. The objective here will be to vary and hopefully enhance the uv 
emission spectrum from the auxiliary discharge, thereby influencing the 
level of photoionlzatlon in the main discharge region. 
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2.0 lASER MODELING 


2.1 Probabilities for V-V Exchange In CO-N ^ and 

2.1.1 CO-N 2 V-V Probabilities 

Calculations of vlbratlonal^lbratlonal energy exchange probabilities 
In pure CO are usually based on the assumption of Independent contributions 
from short range and long range forces. The long range contribution Is 
due to multlpole-multlpole Interactions [1], of which the d Ipole-d Ipole 
term Is dominant. The short range contribution arises from coulomblc 
repulsive forces. 

For molecular systems involving which has no permanent dipole 
moment* long range contributions are usually neglected in calculating V-V 
exchange probabilities. In models of vibrational kinetic processes In 
CO lasers* the CO-N 2 V-V probabilities are typically calculated with 
the short-range expression used for the pure CO calculation [2,3] with some 
modification to obtain agreement with experimental data. An expression of 
this type* based on the Schwartz-Slawsky-Herzfeld theory, was discussed in 
an earlier progress report.* Attempts to use this expression in calcu- 
lating CO-N 2 probabilities* however* led to poor agreement with experi- 
mental results in either the vibrational quantum number dependence or the 
temperature dependence of the V-V probabilities. Much better agreement 
was obtained by abandoning the SSH theory and using in its place an 
expression obtained by an extension of the theory of H. K. Shin [A]. The 


*Super sonic CO Electric Discharge Lasers, NASA Grant NSG-2050* Semi-annual 
Progress Report for the period November 1* 1974 - April 30, 1975. 
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remainder of this section Is a brief discussion of Shin’s theory and its 
extension, vlth comparisons with experimental data. 

The theory is a semlclasslcal treatment, assuming harmonic oscillator 
vave functions and a three-dimensional classical collision trajectory. 

The overall Interaction potential is assumed to be the sum of four atom- 
atom interaction potentials, which are represented by Morse functions. 

The interaction thus Includes both attractive and repulsive terms. The 
molecules are assumed to be rotationless. A result is obtained only for 
the vibrational transitions involving the ground level and the first 
excited level, l.e., in the case of N 2 and CO, for the processes 

N2(v= 1) + C0(v=0) H2(v=0) + CO(v*l) . 


For the forward process. Shin's result may be expressed in the form 


p(C^N2) ^ 
0,1 
1.0 


.1/2 


(1 79xl0~^)f-^^ expl-^ + ^ 

U./yxxu exp^j^ + 


1/2 


16D 

3it\t 


g 



In this expression (o = - u)__ is the difference in vibrational 

N2 CO 

frequencies of the two molecules, D is the well depth of the atom-atom 
interaction potentials, and g is a factor which arises in the particular 
formulation of the overall interaction potential. For the system of N 2 
and CO, Shin's value for g is 0.712. The qiiantlty x ^ equation (1) 
is related to the difference in vibrational frequencies, as 


X 



TTooa XtI 


2/3 


( 2 ) 


where y is the reduced mass for the collision system and a is the range 
parameter in the Interaction potentials. Shin's expression for the transi- 
tion probability, equation (1), yields values which are in good agreement 
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with experimental data at both low and high temperatures [5]. 

Models of CO vibrational kinetics require a method of estimating the 
probabilities for single-quantum vibrational energy exchange involving 
any vibrational level. A simple means of extending Shin's result to the 
general V-V exchange reaction. 


C0(v“r-1) + N 2 (v=s) X CO(v=r) + N2(v=s-1) , 


is to assume harmonic oscillator scaling of the V-V probabilities, l.e.. 


^r-l,r “ ^0,1 


(3) 


s,s-l 


1.0 


Equation (3) which follows from an assumption that the interaction 

Hamiltonian is independent of the initial states, is exact to second order 

li. . 

1.0 

A further modification is to introduce anharmonlcity effects by 
replacing iia> in equation (1) with the actual energy nonresonance, |AEj , 
and to correct the vibrational frequencies with anharmonlcity factors: 


0 ). 






(4) 

(5) 


These substitutions are ^ hoc , and they are justified only by the better 
agreement with experimental data which is obtained. 


*Note, however, that if this substitution is made in the pre-exponential 
term in equation (1), the probability for an exactly resonant V-V r^xchange 
would be zero. Although there :e no exact resonances for purely vibra- 
tional transitions in CO and N 2 » this perhaps academic difficulty is avoided 
by not making the substitution in the pre-exponential term. 
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A comparison of the probabilities calculated by the above method 
with the room-temperature data of Hancock and Smith [7] shows good agree- 
ment, although the calculated probabilities are generally larger by about 
a factor of two. Better agreement is obtained if the constant in the pre- 
exponential is chosen to fit the room-temperature data for . 

i.’o 

The final expression for the probability for the single-quantum V-V reaction 
C0(v=r-1) + N 2 (v=s) CO(v=r) + H2(v=s-1) 


r 

s 




is then 

^(C^N2> ^ (3.52x10" 
r-l,r 

s,s-l 

The function G(x;T) in this expression is defined as 

G(XiT) - (^) . exp i[(-3x + 8 + ^ 8^)] 


where 


( 6 ) 


(7) 


X = 





( 8 ) 


and the energy nonresonance is given by 

.8 - {.y . - (Ef °> - . (9) 

Equation (6) is the expression used to calculate the probabilities for CO-N 2 
V-V exchange in C0ED2. Following Shin, the values for the potential para- 
meters used in these calculations are 
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a - 0.2 

D ■ 1.385 X lO"^^ ergs . 

Figures 1 and 2 are a comparison of the available low- temperature data 

[6-10] with the probabilities calculated from equation (6). Also plotted 

in these figures are the probabilities which would be obtained from the 

short-range SSH expression, with the constant In the pre-exponent lal 

chosen to match the room-temperature data for ^2^ . 

0,1 

1.0 


2.1.2 N 2 -N^ V-V Probabilities 

A similar approach is taken to calculate the probabilities for V-V 
exchange in pure N 2 . Shin’s approach leads to a value of zero for exact 
resonance, so that in this case the process 

N2(v=0) + N2(v=2) 2N2(v*1) 


is used as a basic process to which the vibrational quantum number scaling 
is applied. The calculation of the probability for this reference V-V 
process is the same as Shin's development for CO and although an addi- 
tional factor of •/! ap^cirs in the evaluation of the scattering matrix 
element due to the choice of Initial states. The extension to the general 
single-quantum V-V process, 

N2(v=r-1) + N2<v=s) -»• N2(v=r) + N2(v=s-1) , 


is made in exactly the same fashion as for the CO-N 2 probabilities. 
The final result for the N 2 -N 2 V-V probability is expressed by 


(N2-N2) 

r-1, 


r 


s,s-l 


(3.47x10"^) 




( 10 ) 
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Figure 1. Room-temperature probabilities for the CO-N 2 V-V process 
C0(v) + N2 Cv* 0) -»■ CO(v-l) + N2 (v» 1). The calculated 
probabilities are fit to a value of 1.9 x 10"^ for 

oil 
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Figure 2. Comparison of the calculated probabilities With the low- 

temperature experimental results of Stephenson and Mcsburg 
for the V-V process. CO(v*l) + N2(v*0) C0(v»0) + N 2 (v«l). 
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The definitions of the terms in this equation are the same as in 
eqxiatlon (6) for the CO-N 2 case. The constant g has a value of 0.731 
for N 2 -’N 2 » and the values used for the potential parameters are 

a - 0.2 

D • 1.263 X 10 ergs . 

Calculated probabilities for the V-V excitation of ground-level M 2 

are plotted in Figure 3. The value calculated by Berend and Benson [11] 

for p^^2”^2) at X » 300 K is also shown on this plot Their calcu- 

2.*1 

lation is based on a rigorous two-dimensional classical model, which 

« 

Includes the effects of molecular rotation and anharmonlcity. 

2 . 2 Calculations of CO/Mq Laser Kinetics Using C0ED2 

With the inclusion of N 2 in the model of CO vibrational kinetics, 
C0ED2 has become a useful tool In studying the relative Importance of 
direct electron pumping and CO-CO and CO-N 2 V-V transfer processes in 
shaping the CO vibrational distribution in fast-flow continuous electric 
discharge systems. Several calrclations have been made for mixtv'*es of 
N 2 and CO to assess the effects of the relative N 2 /CO concentrations on 
the CO vibrational kinetics and small-signal gain. These calculations 
assume that the discharge parameters, n^ and T^ , are not affected by 
changes in the relative concentrations of the two molecular species and 
thus Isolate, perhaps artificially, the vibrational kinetics from the 
overall discharge kinetics. Such predictions are useful, however, in that 
they can help Identify phenomena in the vibrational energy transfer pro- 
cesses which could be used to help meet particular objectives in a laser 
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(Wi-Wt) 



design, e.g. high gain in certain vibrational bands, changes in the delay 
tine to naxinun gain, etc. 

Calculations for two mixtures of CO and N 2 are discussed in this 
section. The geometry and flow conditions for these calculations corres- 
pond to the small supersonic double-discharge system at Stanford. The 
flow channel in this system has a two-dimensional nozzle with an area ratio 
of 17.4. The total length of the test section is 30 cm, including the 
20- cm discharge section and a lO-cm section downstream of the discharge. 

The sidewalls of the test section diverge uniformly to allow for boundary 
layer growth, with entrance and exit widths of 1.9 cm and 2.86 cm, respec- 
tively, The top and bottom walls of the discharge channel are formed by 
the flushed-mounted main discharge electrodes, which have a constant 
separation of 1.2 cm. Continuous supersonic flow in the test section is 
maintained by a diffuser and pump system with a volumetric capacity of 94 
llters/sec. 

For the calculations reported here stagnation conditions of 300 K 
and 4.06 atm were assumed. The calculated temperature and pressure in the 
test section for these conditions are approximately 65 K and 14 torr. The 
pressure calculation, which reflects boundary layer corrections to the one- 
dimensional isentropic flow equations, compares favorably with measured 
static pressures in the test section. Gas heating effects due to the 
discharge are included in the calculations. 

Discharge conditions were assumed on the basis of an average current 

2 "16 2 
density of 10 mA/cm^ and an estimated E/N of 2 x 10 V • cm for the 

main discharge. Solutions of the electron Boltzmann equation by Morgan 

and Fisher [12] and Nighan [13] were used to estimate the electron drift 
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velocity and average electron energy. The discharge parameters used In 

10 -3 

the calculations, based on these considerations, were n^ * 2 x 10 cm 

and » 7000 K. If nearly all of the discharge energy is assumed to go 

Into vibrational excitation, these values of electron number density and 

_2 

temperature lead to an energy loading ^ 10 eV/molecule. This energy 
loading is consistent with the discharge performance attained in initial 
supersonic experiments with . ■ 

Stagnation and dlscliarge conaitions for the two calculations are 
sunmiarlzed in Table 1. 

Table 1. Discharge and stagnation conditions 
used in the calculations. 


^0 ‘ 

stagnation pressure 

4.06 atm 


stagnation temperature 

300 K 

% » 

electron number density 

2 X 10 cm 

^e » 

electron temperature 

7000 K 


2.2.1 Case 1, N2:C0 = 90:10 

The first calculation is for a mixture of 90Z N 2 , lOZ CO. Results 
for the CO vibrational distribution and the maximum line-center gain 
coefficients are plotted in Figures 4 and 5. Distance on these graphs is 
defined from the beginning of the 20-cm discharge section. A temperature 
increase in the flow direction of approximately 3K is predicted due to gas 

heating effects. The ca3culated vibrational energy densities at the end 

-2 *3 

of the discharge region are 2.5 x 10 eV/raolecule for CO and 5.4 x 10 

eV/molecule for N 2 • 
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The variation of the vibrational distribution and maxinum gain 
coefficients with distance along the flow indicate that the CO-CO V-V 
processes dominate the vibrational energy transfer for this calculation. 
These effects are seen in the Increase in gain of the higher -lying laser 
transitions at the expense of the gain on the lower transitions, repre- 
senting an upward flow of vibrational energy. The net vibrational energy 
transfer in the inter species V-V processes is from Nj to CO, due to the 
smaller spacing of the CO vibrational energy levels. Most of this trans- 
ferred energy is pumped into the low-lying CO levels, since the energy 
spacings for the lower levels are more nearly resonant with the energy 
spacings. This process is slow due to the small excited -level N 2 popul- 
ations, and the transferred energy is quickly redistributed in CO by the 
faster CO-CO V-V processes. Although the probabilities for near-resonant 
V-V transfer in CO-CO and CO-N 2 are of similar magnitudes, the rates of 
transfer are also dependent on populations of vibrational levels which 
participate in the near-resonant transitions. For CO-CO collisions, these 
levels are adjacent and are rapidly populated through a "ladder "-climbing 
process, whereas for CO-N 2 collisions the near-resonant exchanges involve 
weakly-populated high- lying levels in N 2 . 

Total inversions are predicted at Intermediate locations in the dis- 
charge for the A -► 3 and 5 -*■ A transitions. These inversions lead to 
significant gain for R-branch transitions. The P-branch and R-branch gain 
for the A 3 vibrational band is plotted in Figure 6, for a distance of 
5.5 cm into the discharge. 


15 


Figure 6. P-*branch and R-branch gain coefficients at line center for 
the V “ 4 -► V ■ 3 vibrational band. The distance from the 
beginning of the discharge Is 5.3 cm. 




2.2.2 Case 2, N^rCO - 98:2 

The second calculation Is for a mixture of 98Z and 2% CO. The 
change in CO partial pressure from Case 1 causes a reduction by a factor 
of 5 in the CO-CO V-V transfer rate. This effect is evident in the 
calculated CO vibrational distribution. Figure 7, where the increase in 
population for the highest levels is slower than in the first calculation. 
The most striking difference in the calculated distributions for the two 
cases, however, is in the energy content, or average vibrational energy 

per CO molecule. At the end of the discharge, the vibrational energy . 

-2 -3 

densities for Case 2 are 5.8 x 10 eV /molecule for CO and 6.3 x 10 

eV/molecule for N 2 . in the lO-cm section downstream of the discharge, 

_2 

the CO vibrational energy increases to 7.7 x 10 eV/molecule due to V-V 
transfer from N 2 . 

The higher vibrational energies in Case 2 are due to a greater energy 
storage capacity of the gas mixture. V-V transfer out of N 2 in CO-N 2 
collisions, per molecule of N 2 , is decreased due to the factor of 5 reduc- 
tion in the CO concentration. At the same time, the rate of V-V energy 
transfer from N 2 into CO, per molecule of CO, is actually Increased due to 
the slight Increase in N 2 concentration. The result of these two effects 
is a higher CO vibrational energy, on a per -molecule basis. 

The maximum gain curves for Case 2 are plotted in Figure 8. The over- 
all magnitudes of the small-signal gain are generally lower than in Case 1, 
due primarily to the decreased CO concentration which results in lower 
population difference densities. The high gain for low-lying transitions 
in the early stages of the calculation is due to direct electron pumping. 
The increase in gain for these transitions farther downstream, particularly 
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beyond the discharge region* reflects the energy transfer from N2 into the 
low- lying levels. This effect is of potential Interest* because it may 
produce high gain for the lov-lying transitions in a region downstream of 
the discharge* i^ere the optical quality of the medium is perhaps Improved. 
The phenomenon of delayed low-level gain enhancement due to V-V pumping by 
N2 has also been observed in the CO laser kinetics calculations of Fisher 
[ 14 ]. 
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3.0 LABORATORY STUDIES OF DOUBLE DISCHARGE 

Experimental vork during this six month period has continued to focus 
on the concept of the CW double discharge. In order to better conduct 
tests aimed at elucidating the discharge mechanism and Increasing the 
electrical power loading* we have constructed a second-generation experi- 
mental facility. Although our available pump capacity limits the super- 
sonic channel to the same cross-sectional dimensions employed in our first 
facility, the new facility has several improved design features and much 
greater flexibility. The basic design was described in our previous 
progress report. During the past six months ^.he system has been assembled 
and tested for gasdynamic and electrical performance. 

In summary, our experimental findings are as follows: 

1. The cathode material was found to be an important factor in 
discharge performance, with copper performing better than 
aluminum. 

2. Water cooling of the electrodes and Improved mounting of pins 
has enabled operation at considerably higher pin currents 
than possible with the original test section. 

3. The maxlrnimi electrical power loading (before arcing) was found 
to scale with the total pin current. At pre' ent we are power 
supply limited and cannot determine the xlmum power loading 
achievable with total pin currents larger than 300 mA. We 
regard this effort to determine the dependence of maximtuB power 
loading with pin current to be a critical element of our study, 
and we are hopeful that the power loadings required for lasing 
can be achieved by operating with increased pin currents. 
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A. Prellnlnary teats with individually variable ballast resistors 
in the pin circuitry have confirmed that control of the indi- 
vidual pin currents can be used to improve discharge imlformlty. 
Discharge uniformity and electrical power loading are of course 
both important in a practical device. 

Recently we have ordered a new power supply to extend our operating 
range to higher power loadings. In>the meantime, ve have also made a 
major modification to the test section to allow injection of foreign 
gases into the cathode and anode boundary layers. This modification is 
motivated by our recognition of the role of photoionization in stabilizing 
the main discharge. Varying the gas in the vicinity of the cathode, where 
the photoionlzing photons are created, should enable variation of the 
emission spectrim and hence of the volume ionization level in the core 
flow. At present, we are initiating tests to check the gasdynamlc per- 
formance with injection. Subsequently, when the new power supply is 
installed, electrical tests will be conducted. 
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4.0 RESEARCH FOR NO\fEMBER 1, 1975, TO APRIL 30, 1976. 


During the forthcoming six-month period we plan to work on the 
following problems: 

A* Laser Modeling 

1. Parameter studies of system performance as a function of 

key variables such as gas composition, temperature, pressure, 

power loading and kinetics coefficients. 

2. Addition of an option to Include an optical cavity, if needed. 

3. Development of a physical model for double discharge. 

B. Laboratory Studies 

1. Work to improve system performance, particularly the discharge 

uniformity and the power loading of the gas. Ideas Include: 

a. Individual variable ballasting of pins in auxiliary 
discharge. 

b. Segmentation of anode for axial control of E/N in the 
main discharge. 

c. Studies at Increased pin and main discharge currents, 
made possible with a new power supply. 

d. Electrode modifications, including rounding downstream 
edges of cathode and anode and masking anode sides to 
reduce breakdown through the warm sidewall boundary 
layers. 

e. Injection of gases into the cathode boundary layer, 
to stimulate UV emission or provide better heat dis- 
sipation. 

f . Variation in number, location and s^ze of pins. 
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2. ^rk to loprove understanding of double discharge: 

a. Investigate the role of photoionization 

b. Investigate the role of electrode naterials and their 
characteristics such as work function* cathode voltage 
drop and naterial cleanliness and purity. 

c. Develop physical model of double discharge. 

During the forthcoming year* we anticipate that our NASA- sponsored 
research will benefit from funding recently received from the National 
Science Foundation. These funds will allow us to Improve and enlarge our 
experimental facilities for stiidies of discharges in supersonic flows. 

The overlap between the NASA- and NSF-sponsored research should benefit 
both programs* and in particular the Increased size of our next test 
sectlcm* once Installed* will enable more practical experiments than 
presently possible. 
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